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Abstract 



In this paper, we establish two new estimates concerning positive solutions of a linear 
parabolic equation. More precisely, we give explicit estimates of the gradient of logarithm 
fV^ of a positive solution via the uniform bound of the logarithm of the solution. Moreover, we give 

PLh 1 a generalized version of Li-Yau's estimate. Our proof is based on the link between PDE and 

^ (-h quadratic BSDE. 

> 

1 Introduction 

In this article, we study positive solutions u of a linear parabolic equation 
J> ^L-^ju = in(0,oc)xM, (1.1) 

in, 

where M is either the Euclidean space M. n and L is an elliptic or sub-elliptic operator of second- 
order L = I Y^=i Aa + Aq, {Aq, • • • , A m } is a family of vector fields on W 1 , or M is a complete 
manifold of dimension n with Riemannian metric (gij), and 2L is the Laplace-Beltrami operator 

A ^ ^ l~ l 3 ^ 

a; 

where 5 denotes the determinate of (<7ij) and (g 1 - 7 ) is the inverse of the matrix (gij). 

The well-posedness and the regularity theory for (11.11) are parts of the classical theory in partial 
differential equations, see [18], [13J and [20J for details. On the other hand, it remains an interesting 
question to devise precise estimates of a solution u in terms of the (geometric) structures of (jl.ip . 
There is already a large number of papers devoted to this question. Among many interesting 
results, let us cite two of them which are most relevant to the present paper. The first result is a 
classical result under the name of semigroup domination, first discovered by Donnelly and Li [8], 
which says that if the Ricci curvature is bounded from below by C, then 

|VP t u | < e~ Ct Pt\Vu \ for u € C£(M) (1.2) 
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where (Pt)t>o is the heat semigroup on M, so that the left-hand side is the norm of the gradient of 
u(t, •) = PtUQ a solution to the heat equation 

1 d \ 

-A-— Ju = in(0,oo)xM (1.3) 

with initial data u(0, •) = no, while the right-hand side Pt|Vuo| is a solution of (|1.3p with initial 
data |V«o|- The second result is Li-Yau's estimate first established in [19] . If the Ricci curvature 
is non-negative, and if u is a positive solution of (jl.3p then 



IV log ul 2 - 2— log u < - fori>0. (1.4) 
dt t y ' 

In fact, in the same paper [19J, Li and Yau also obtained a gradient estimate for positive solutions in 
terms of the dimension and a lower bound (which may be negative) of the Ricci curvature, though 
less precise. Their estimates in negative case have been improved over the years, see for example 
[25], [26], [2] and pQ. 

In this paper we prove several new gradient estimates for the positive solutions of (jl.lj) . 

Theorem 1.1 Let M be a complete manifold with non-negative Ricci curvature. Suppose u is a 
positive solution of U.3\) with initial data uq > 0, then 

|Vlogn(t,-)| 2 < ^Hloguolloo fort>0 (1.5) 
where \ \ ■ ||oo denotes the L°° norm on M. 

Indeed we will establish a similar estimate for the heat equation with a sub-elliptic operator, 
under similar curvature conditions, and indeed we will establish a gradient estimate for a complete 
manifold whose Ricci curvature is bounded from below. We should mention that the kind of 
estimates such as (|1.5|) is known in the literature, and it is known that 

C 

I Vit(i, -)| < — p| \uq\ |oo for t > 
yt 

for some constant C > which depends on the manifold M in a complicated way. The interesting 
feature in our estimate (|1.5p is that there is no unknown constant, and (II. 5h does not involve the 
dimension. 

Theorem 1.2 Let M be a complete manifold of dimension n with non-negative Ricci curvature. 
Suppose u is non-negative solution to the heat equation of U.3\) with initial data uq > 0. 
1) If C 6 [0, oo] such that — Aloguo < C, then 



8. ^ C 

.— logu < -. 

dt ~ i-C + 1 



iVlogul 2 - 2— logti < -. fort > 0. 



2)IfCe (0,oo) such that -Alogu > -C, then 



d C n 

|Vlogn| 2 -2— logu> -— -j- forte[0,-). 

n 
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By setting C = oo in 1) we recover Li-Yau's estimate (|1.4p . 

The novelty of the present paper is not so much about the new gradient estimates in Theorem 
ll.ll and Theorem ll.2l though we believe that they are not known previously even in the case that the 
underlying manifold is Euclidean space, what is interesting of the present work is the approach we 
are going to develop in order to discover and prove these gradient estimates. Our approach brings 
together with the martingale analysis to the study of a class of non-linear PDEs with quadratic 
growth. Of course the connection between the harmonic analysis, potential theory and martingales 
is not new, which indeed has a long tradition, standard books may be mentioned in this aspect, 
such as [9], [10], [IT] and etc., what is new in our study is an interesting connection between the 
BMO martingales and positive solutions of the heat equation (jl.3p . 

To take into account of the positivity, it is better to consider the Hopf transformation of a 
positive solution u to (|1.3p . i.e. / = log it, then / itself solves a parabolic equation with quadratic 
non-linear term, namely 

1 . d\ „ 1,_„, 2 



2 A--l/ = --|V/|^ in[0,oo)xM. (1.6) 

The preceding equation (|1.6p is an archetypical example of a kind of semi-linear parabolic equa- 
tions with quadratic growth which has attracted much attention recently associated with backward 
stochastic differential equations, for example Kobylanski |16| . Briand-Hu [6], Delbaen et al. [7] and 
etc. 

The main idea may be described as the following. Suppose / is a smooth solution of the non- 
linear equation (|1.6p . and Xt = Bt + x where B is a standard Brownian motion on a complete 
probability space. Let Y t = f(T - t,X t ) and Z t = (Z l t ) where Z\ = V*/(T - t,X t ), V 4 is the 
covariant derivative written in a local orthonormal coordinate system. Then, Ito's lemma applying 
to / and X may be written as 

Y T-Yt = Y.J t KdBl \Z s \ 2 ds. (1.7) 

On the other hand, it was a remarkable discovery by Bismut [lj (for a special linear case) and 
Pardoux-Peng [23] that given the terminal random variable Yt S L 2 (VL, Tt, f), there is actually a 
unique pair (Y, Z) where Y is a continuous semimartingale and Z is a predictable process which 
satisfies f)1.7|) . The actual knowledge that Z is the gradient of Y may be restored if Yt = fo(Xx)- 
The backward stochastic differential equation (jl.7jl with a bounded random terminal Yt, which has 
a non-linear term of quadratic growth and thus is not covered by Pardoux-Peng [23] , was resolved 
by Kobylanski [16] . Observe that the martingale part of Y is the Ito integral of Z against Brownian 
motion B (which is denoted by Z.B). It can be shown that, if Y is bounded, then Z.B is a BMO 
martingale up to time T, so that the exponential martingale 



£(h(Z).B) t = exp 



n r-t If* 

/ h\Z s )dB\-- / \h(Z s )\ 2 ds 
~1 Jo 1 Jo 



is a uniformly integrable martingale (up to time T), as long as h is global Lipschitz continuous. The 
main technical step in our approach is that, due to the special feature of our non-linear term in (jl.6p . 
we can choose h l (z) = z l (one has to go through the detailed computations below to see why this 
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choice of K l is a good one), and making change of probability measure to Q by = £(1i(Z).B)t, 
then, under Q, not only Z.B is again a BMO martingale (where B is the martingale part of B 
under the new probability Q), but also t — > \Zt\ 2 is a non- negative submartingale. Next by utilizing 
the BSDE (jl.7j) . we can see the BMO norm of Z.B under Q is dominated at most 2-^/1 |y| |oo, that 
is 



E®1 I \ZA 2 ds 



T t > <4||y| 



Finally the sub- martingale property of \Z t \ 2 allows to move \Z S \ 2 (for s E (t,T)) out from the time 
integral on the left-hand side of the previous inequality, which in turn yields the gradient estimate. 

Let us now give a heuristic probabilistic proof to Theorem 11.21 to explain from where such 
estimates come from. Let / = log-u, and G = —A/. Then one can show that 

G = |V/| 2 -2/ i 

(where ft stands for the time derivative J^/ for simplicity). Moreover, G satisfies 

(L-^)G = ±G 2 + H, 
ot n 

where 

H = (| VV/| 2 - -G 2 ) + 2Ric(V/, V/), 
n 

and H > 0. We suppose here G > 0. Consider the BSDE: 

dY t = Z t dB t + -Y 2 dt, Y T = G(0, x + B T ). 
n 

Then 

Y t >G(T-t,x + B t ). 

Setting 



i - 1 v- Zt 

then (U, V) satisfies the following quadratic BSDE: 



1 \v t \ 2 
dU t = — dt + v t dB t + 

n U t 



Using BMO martingale techniques, one can prove that there exists a new probability measure 
under which Bt = Bt + J Q * ^f-ds is a Brownian motion. Hence 

dU t = --dt + V t dB t , 
n 

from which we deduce that U = | + E Q [£/" T ], and 

1 



Yr 







i 
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which yields the estimate in Theorem 11.21 

Even though the above heuristic proof is probabilistic (which can be made rigorous), we prefer 
to give a pure analytic proof in the last section. 

The paper is organized as follows. Next section is devoted to some basic facts about quadratic 
BSDEs including BMO martingales. Section 3 establishes the gradient estimates for some lin- 
ear parabolic PDEs on Euclidean space, while Section 4 establishes these estimates on complete 
manifold. Last section is devoted to establish a generalized Li-Yau estimate via analytic tool. 

2 BSDE and BMO martingales 

Let us begin with an interesting result about BSDEs with quadratic growth. The kind of BSDEs 
we will deal with in this paper has the following form 

m m 

dY = Y J Z j F j {Y,Z)dt + s ^Z j dB j , Y T = i, (2.1) 

i=i j=i 

with terminal value £ G L°°(Q, J-j>, IP) which is given, where B = (B ,B m ) is a standard 
Brownian motion, (J r t)t>o is the Brownian filtration associated with B, and F 3 are continuous 
function on R x R m with at most linear growth: there is a constant C\ > such that 

\F(y, z)\ < Ci(l + \y\ + \z\) V(y, z) G R x R m . 

According to Peng [24] and as we have seen in the Introduction, if u is a bounded smooth 
solution to the following non-linear parabolic equation 

^ u + J2Fi(u,Vu)^ r = ^Au in[0,oo)xR m (2.2) 

3=1 3 

with initial data Uo, then Yt = u(T — t,Bt + •) and Zt = Vu(T — t,Bf + •) is a solution pair of (|2.ip 
with terminal value Yt = uq{Bt + •)• The special feature of (|2.2|) is that the maximum principle 
applies, which implies that global solutions (here global means for large t) exist for the initial value 
problem of the system as long as the initial data is bounded (though, this constraint can be relaxed 
a bit, but for the simplicity we content ourself to the bounded initial data problem). The maximum 
principle implies that as long as u is a solution to (12. 2p then |?x(x,t)| < ||uo||oo- Therefore, if the 
initial data uq is bounded, and F J are global Lipschitz, then, according to Theorem 6.1 on page 
592, [H], u exists for all time, and both u and Vit are bounded on R m x [0,T]. 

The maximum principle for (|2.1|) however remains true even for a bounded random terminal 
value (so called non Markovian case), which in turn yields that the martingale part of Y is a BMO 
martingale. This is the context of the following 

Proposition 2.1 Suppose that £ G L°°(Q, Tt, P)- There exists a unique solution (Y, Z) to h2. 1\) 
such that Y is bounded and M = Z.B is a square integrable martingale. Moreover M = Z.B is a 
BMO martingale up to time T, and 

im*)iioo < iieiioo vtG[o,r]. 
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Proof. The existence and uniqueness is already given in [16]. The fact that M = Z.B is a 
BMO martingale up to time T is proved in [22]. Then there exists a constant C2 > such that 

r rT 



E 



\ZJ ds 



< Co 



Let N t = Ej=i Jo FJ (Ys, Z s )dBl Since 
(N,N) T -(N,N) t 



Y,\F'(Y s ,Z s )\ 2 d S 

j 

< J Cf(l + \Y S \ + \Z s \) 2 ds, 



so that there exists a constant C3 > such that 

E{{N,N) T - (N,N) t \Tt} <C 3 . 

Therefore N is a BMO martingale. Hence the stochastic exponential £(—N) is a martingale up to 
T. 

Define a probability measure Q on (£7, Tt) by dQ/d¥ = £(—N)t- Then, according to Girsanov's 
theorem Bt = Bt + (iV, B)t is a Brownian motion up to time T under Q, and (Y, Z) is a solution 
to the simple BSDE 

dY t = Z t .dB t 



under the probability Q, whose solution is given by 



Y t = E^{£|JU = K{£(-N) T £(-N)-^\T t } for t < T. 



(2.3) 



It particularly implies that ||lt||oo < HClloo- 



3 Stochastic flows and gradient estimates 

Let Aq, A\, • • • , A m bem + 1 smooth vector fields on Euclidean space R n , where n is a non-negative 
integer. Then, we may form a sub-elliptic differential operator of second order in M. n : 



1 m 

L=-Y,A 2 a + A (3.1) 



2 

a=l 

here we add a factor \ in order to save the constant v2 in front of Brownian motion which will 
appear frequently in computations in the remaining of the paper. Our goal is to devise an explicit 
gradient estimate for a (smooth) positive solution u of the heat equation 

L- — \u = 0, on (0, 00) x R n , (3.2) 

by utilizing the BSDE associated with the Hopf transformation / = logn, which satisfies the 
semi-linear parabolic equation 

m 

on (0, 00) x R n . (3.3) 



a=l 
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3.1 Stochastic flow 



The first ingredient in our approach is the theory of stochastic flows defined by the following 
stochastic differential equation 

m 

dip = Ao(tp)dt + ^ A a(p) ° dw a , ip(0, •) = x (3.4) 

a=l 

where od denotes the Stratonovich differential, developed by Baxendale [3], Bismut [5], Eells and 
Elworthy [12j . Malliavin [21], Kunita [T7j and etc. The reader may refer to Ikeda and Watanabe [13] 
for a definite account. To ensure the global existence of a stochastic flow, we require the following 
condition to be satisfied. 

Condition 3.1 Let A a = X}j=i ^"afj • Assume that A 3 a have bounded derivatives. 

By writing (|3,4p in terms of Ito's stochastic integrals, namely 

m 

(p)dt + AUv)dw a , <p(0, •) = x, (3.5) 

a=l 

where (and thereafter) Einstein's summation convention has been used: repeated indices such as 
/ is summed up from 1 up to n. The existence and uniqueness of a strong solution follow directly 
from the standard result in Ito's theory, which in turn determines a diffusion process in M. n with 
the infinitesimal generator L. 

In fact, more can be said about the unique strong solution, and important consequences are 
collected here which will be used later on. Suppose w = (wt) is a standard Brownian motion (started 
at 0) with its Brownian filtration {J~t)t>o on the classical Wiener space (f2, J 7 , P) of dimension m, 
so that w = (wt)t>o is the coordinate process on the space f2 of continuous paths in M. m with initial 
zero. Then, there is a measurable mapping ip : R + x 17 x K n — > W 1 and a probability null set jV, 
which possess the following properties. 

1. w — > (p(t, w, x) is JVmeasurable for t > and i£K™, and (p(0, w, x) = x for every w £ $7\W 
and i £ I", 

2. t — > (p(t,w,x) is continuous, that is ip(-,w,x) G C(M + ,lR ra ), for w G $7 \ TV and x G M. n . 
t —> (p(t, - ,x) is a continuous semimartingale for any x G M n . 

3. x — >■ cp(t, w, x) is a diffeomorphism of M. n for each w G f2\W and t > 0. That is x — )• </?(i, u;, x) 
is smooth and its inverse exists, and the inverse is also smooth. 

4. The family {<p(t, -,x) : t > 0, x G M n } is a stochastic flow: 

9?(t + s, u>, x) = <^(t, 9 s w, cp(s, w, x)) 

for all t, s > 0, x G M n and m)£H \ A/", where ^ s : O — > Q is the shift operator sending a path 
to to a path 9 s w(t) = w(t + s) for t > 0. 

5. For each x G M n , <£>(£) = (p(t,-,x) (or denoted by tp(t,x)) is the unique strong solution of 



g4 

9x J 
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6. Let JUt, w, x) 



dip 1 (t ,w ,x) 



-Q-j — for i,j < n. Then Jj(0, w, x) = 5) and J solves the following SDE 

(3.6) 



ft A 1 J—* ft A 1 

Oi=l 



and its inverse matrix K = J = (K 1 -) solves 



dK) 



ftA 1 ftA 1 
a=l 



(3.7) 



In our computations below, we have to use Ito's integrals rather than Stratonovich's ones. 
Therefore we would like to rewrite (|3.6l I3.7P in terms of Ito's differential, so that 



dJ) 



a=l 



and 



dK) 



dx l 

a=l x 



^ ftA 1 

a=l 



a _|_ dA k a dA l a \ 
dx l dx k dx l dx k J 



{<P)J l jdt 



dxi 
8A l n 1 



dxi 



a=l 



. k d 2 A l a dA k dA l Q 



' dxWx k dxi Qx k 



(tp)dt. 



(3i 



(3.9) 



3.2 Structure assumptions 

We introduce some technical assumptions on the structure of the Lie algebra generated by the 
family of vector fields {A$,Ai,--- ,A m }, in addition to Condition 13.11 Recall that A a = Aa-?—, 

and A J R , A 3 R etc. are the corresponding coefficients in Lie brackets 



so that 



d 

[A a ,A p ] = A 3 afi — p [A a , [Ap,A 



Al a = A 



j _d_ 

A a,f3,V dx j etC '' 



dAi . QAl 



$ a dx i 



A\ 



dx % ' 



etc. Let 



A k 



AiA 



d 2 A k 



AiAt 



d 2 A k 



/9 dx i dxi 13 a d&dxi 

. dA i dA k . dAi dA k 
-2Ai f a „ f +AI P 13 



dA k dM 



dxi Q x % 



13 dxi dx % 



dx l dx'i 



R k a = J2 A h,« = T,i A ^ A e> A ^ k ■ 

0=1 /3=1 



(3.10) 
(3.11) 
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Condition 3.2 There is a constant C\ > such that for any £ = (£«)i< n , Op = (Oi t p)i< n S 
(ft = 1, ■ ■ ■ , mj, ZioZds i/iat 

m / n \ ^ m / n \ / n \ 

E E4*U +2 E E^4,J E^M 

a,/3=l \fe=l / a,/8=l \fc=l / \i=l / 

m / n \ / n \ 



a,/3=l \fe=l / \i=l 

m / n \ ^ 



a=l \fe=l 



?.e. 



a,/3=l 

m 

> -dE^-o 2 - 

a=l 

Condition 3.3 T/iere zs a constant C2 > suc/i i/icti /or any £ = (£i)?<n G M n 

n I m m \ 

E 6 £(4*£ + 2 44«) + E 4,a4,a & 

i,fc=l y«=l a,/3=l / 

m / n \ ^ 



a=l \fc=l / 

The above conditions are satisfied if ^4 is elliptic. Let us now suppose the following Frobenius 
integrability condition: there exist some bounded smooth coefficients c^ Q (x), such that 

m 

A3,a = E c k<*^' ^ = 0,l,...,m, a = l,...,m. 
1=1 

In other words, the Lie brackets .A^q, /3 = 0, 1, . . . , m, a = 1, . . . , m, must lie in the linear span 
of Ai, . . . ,A m . Then the conditions (3.2) and (3.3) are satisfied. 
Indeed, 

r) dA^ 
[Ap,[Ap,A a ]Y = 4—^Aj) - c^Aj-^- 

- J A j 1 dC ^ a Ai A j 
dx l 



This means that [Ap, [Ap,A a ]] also lies in the linear span of A\, . . . ,A m , and the conditions (3.2) 
and (3.3) are easily checked. 
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3.3 The density processes Z a 

Let us consider a smooth solution / to the following non-linear parabolic equation 



£-glJ/ = M/.4>/). onR+xl!" (3.12) 

where h is a C 1 -function on R x R m , though our archetypical example is / = logu and u is a 
positive solution to equation f|3. 3|> . 
By Ito's formula 

Y t = Y T -Y] / Z a dw a - / /i(Y, Z)ds (3.13) 

where 

Y(t,w,x) = f(T -t,tp(t,w,x)), 
Z a (t,w,x) = (A a f)(T-t,ip(t,w,x)) 

for a = 1, • • • , m, and Z = (Z a ). The arguments w and / or x will be suppressed if no confusion 
may arise. Equivalent ly 

m 

= ^ + z)dt. (3.14) 

a=l 

Our aim in this part is to show that Z is an Ito process, and derives stochastic differential 
equations for Z (which in turn gives its Doob-Meyer's decomposition). 

It is clear that both Y and Z a (a = 1, • • • , m) are continuous semimartingales. Taking deriva- 
tives with respect to x l (i = 1, • ■ • , n) in the equation ()3.14j) one obtains 

m / m \ 

dYi = z ?dw a + h y (Y, Z)Yi + h * a (X, z ) Z i dt ( 3 - 15 ) 

a=l \ a=l J 

where 



and 



y = --y and Zf = — Z Q , i = l,---,n 

ox 1 ' ox 1 



^-h(y,z) , h Za = ^- 
oy oz, 

On the other hand, by definition, 

dx* lt " > dif- 

so that 

{t,<p(t,;x))=K l k (t,;x)Y l {t,;x). 



h y = —. h (y' z )' h Za = ^—h(y,z), a = l,---,m. 



Yi(t, -,x) = -K~[Y(t, ; X) = -*y(t, ip(t, ; x))j( (t, -,x) 



dip 

It follows that 



= Afait^xfiKktrrfYfa^x) (3.16) 
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which implies that Z is a continuous semimartingale. The equation (|3.16p is not new, and has been 
used by many authors in different contexts. 

We next would like to write down the stochastic differential equations that Z a must satisfy by 
using the relation (j3. 16|) . which is however an easy exercise on integration by parts. Indeed, we 
have 



f)Ak 

dZ a = Y lK {-^-^W +Y l A k a ^)dK{ + A k a { V )K{dY l 

r)A k r)A k 
+Y l ^(v)d&,Kl)+Kl-^-(<p)d&,Y l ) 

+A k Md(K l k ,Y l ) + ly^^^^jd^V) • 



(3.17) 



Using the SDEs (13.41 13 . 9|) and the BSDE (13.151) . through a length but completely elementary 
computation, we establish the following Doob-Meyer's decomposition for Z 



dZ a = Y^U^Jdw? + h Zp {Y,Z)di) +K l k J2 A Xa Z l dt 



+K i k Y l 



1 m m 
A 0,a + ~ E 4,/3,a + A k a hy(Y, Z)-J2 A«K ( y ' Z ) 



0=1 



dt (3.18) 



(3.19) 



where repeated indices are added up from 1 to ra, 

U a ^ = K l k (A^ a Y l + A k a zP 
and 

A% a = [Ap,A a ] k , A k M>a = [A p , [A p , A a ]] k . 
We are now in a position to work out the Doob-Meyer's decomposition for 

m 

\z\ 2 = Y,\z a \ 2 

a=l 

which simply follows from Ito's formula and (|3,18p . In order to simplify our displayed formula, we 
introduce the following notations: 



3=1 3=1 



for i, k = 1, • • • , n and f3 = 1, • • • , m, so that 



m. I n 



z a = £4& and \ z \ 2 = E 

j=l a=l \j=l 



(3.20) 
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Let 

dw p = dw? + h Zfj (Y, Z)dt 
which is a Brownian motion under probability Q with the Cameron-Martin density 



tflP 





m „ 


= exp 









m ft 1 j-t m 

e/ k^z)^-\ y: 

8=1 J ° 2J ° 8=1 



\h Zg (Y,Z)\ 2 ds 



Then, by an elementary computation, 

m 

d\Z\ 2 = 2 (z% k A% a ^^ a A k J k>p )dw^ 



a,P=l 

m m 

+ E A k a e k ^j hP + 2i k E (A% a Ai + 4 >a A k a )e^ 

a, 8=1 a,0=l 
m m 

+2 h y {Y,Z)Y,A\A k a - Y, A\A k ta h Zp {Y,Z) & k dt 

a=l a,0=l 

m m 

+ Y A i(Ri + 2A 0,a)+ E 4,«4,«\ tekdt. 
a=l a,8=l 

Lemma 3.4 Ifh(Y,Z) = h(Y, \Z\ 2 ), then 

m 

d\z\ 2 = 2^(n t 4+^^)^ 



dt 



a, 8=1 



+ 



+ 



E A k a 9 k ^A i J i]P + 2i k E {4,* A « + 4,« A *I«>.; 

a,8=l «,/3=l 



dt 



a=l 



E 4 (R\ + 1A\ a + 2A k h y (Y, \Z\ 2 )) + E 4,a A 

a,8=l 



k 

8, a 



Proof. In this case 



E &^4,aM y > z )& = 2 E ^^44, a 4^ 

i,k,a,@ i,k,a,8 

= -2 e h%44,«4££k 



i,k,a,/3 



so that 



E SiAiA% a h Z0 (y,z)s k = o 



i.k.a,. 



and thus (^231) follows directly from (13321) . 



(3.21) 



(3.22) 



fc&di. (3.23) 
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3.4 Gradient estimates 

Recall that / is a smooth solution to the non- linear parabolic equation (|3.12|) . where the nonlinear 
term h(Y, Z) has at most quadratic growth. In order to devise explicit estimate for A a f (a = 
1, • • • , m), we assume the following condition to be satisfied. 

Condition 3.5 h(y,z) depends only on (y, \z\ 2 ), i.e. there is a continuously differentiable function 
denoted again by h so that h(y,z) = h(y, \z\ 2 ), and we assume that 

dy 



d\Z\ 2 > -K\Z\ 2 dt + 2[Z a i k A k 0a + i l A l a A k a e k Adwf i (3.24) 



Then, under Conditions 1) - 4), according to (|3.23p . we have 
where K = C\ + C 2 . 

Lemma 3.6 Assume that Conditions i3.ll [gH \EM are satisfied. Then M t = e Kt \Z t \ 2 is 
submartingale under the probability Q (up to terminal time T): 

E Q {M t \F s } >M S VO < s < t < T. (3.25) 

Proof. By Ito's formula 

dM = KMdt + e Kt d\Z\ 2 
so that, for any < s < t < T, we have 



M t 



M s = K f M r dr+ f e Kr d\Z\ 2 

J S J S 

ft m . 

> 2 J e Ks {Z a i k A k ^ + i i ^ a A k a e k ^dw^ 



which yields m 

We are now in a position to prove the following gradient estimate. 

Theorem 3.7 Assume that Conditions L3.ll 1 3. \3.3\) are satisfied. Then 

m AK 

\A a logu(t,x)\ 2 < —K(T—t) II ^ 

^0 1 1 oo 

(3.26) 

a=l 

for any positive solution u of 113. fy) . 



Proof. Apply the computations in the preceding sub-section to / = log u, and h(y, z) 
Then, under the probability Q (defined by (|3.2ip ) 



1 1 „\2 



Y t = Y T - Y, / Z a dw a + / 



J2z a h Za (Y,Z)-h(Y,Z) 



.a=l 



ds, 
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so that 



and therefore 



m T ^ 
Yt = Y T -Y j Z a dw a --\ \Z s \ 2 ds 



E^i^j^\Z s \ 2 ds T t \ = E®{Y T -Y t \T t } 

< 2||y r ||oo<2||loguo||oo 
On the other hand Mt = e Kt \Zt\ 2 is a submartingale, thus one has 

E Q {\Z S \ 2 \ F t ) > e K ^\Z t \ 2 Vs e [t,T] 

so that 

eQ \\f t W ds |^} * \[ z K{t - s) \Z t \ 2 ds 

~ 2K lZtl ■ 

Putting (ET271 [3T28]) together to obtain 



(3.27) 



(3.28) 



\Z t V < 



1 _ e -K(T-t) 



log lto| 



which yields (pT26j) . ■ 

In general, we may proceed with / = ip{u) where tp is a concave function, and u is a positive 
solution to (jl.ip so that / solves (|3.3p with 



h(y,z) 



2 W~\y))\ 2 

We can proceed as above. Under the probability Q 



m ,. T T 

Y t = Y T -Y^ Z a dw a + \ 
a=1 Jt Jt 



J2z a h Za (Y,Z)-h(Y,Z) 



,a=l 



ds, 



so that 



y-y, v r z a dw a + 1 r isrim 

Z-~/ I, Z dw +2j t FW 



a=l 



\Z s \ 2 ds. 



It is important to note that if 



V; (3 V < 2|^T 
then h y (y, z) > 0, so that from Lemma 2.2 in [7j 

fT 



E K 



\Z s \ 2 ds 



^ <4||y T ||^<4|| log^ll 



2 

oo 



(3.29) 
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T 

EM / \Z s \ 2 ds 



Z t \ 2 . (3.30) 



On the other hand M< = e Kt \Zt\ 2 is a submartingale, thus one has 

E Q {\Z S \ 2 \ T t ) >e K ^\Z t \ 2 VsG [t,T 

and therefore 

1 _ e ~K{T-t) 

K 

Putting ([3T291 13301) together to obtain 

\Zt\ < l _ e _ K{T _ t) \\\ogv Ki \\ 00 
which yields the following estimate. 

Theorem 3.8 Assume that Conditions h3. 11 \3.HH \3.3\) are satisfied. Moreover, tp is concave and 
satisfies: 

V> ( V < 2|^"| 2 • 

Then 

m AK 

\A a ^{u{t,x))\ 2 < 1 _ e _ K(T _ t) \Mu )\\ 2 oo (3.31) 

for any positive solution u of 13. 

4 Heat equation on complete manifold 

In this section, we study positive solutions of the heat equation 

1 d \ 

-A - — \ u = in [0, oo) x M (4.1) 

£ CJJ/ J 

where M is a complete manifold of dimension n, A is the Beltrami-Laplace operator. In a local 
coordinate system so that the Riemann metric ds 2 = gijdx % dx^ and 

1 ^ d d 
x fq ^ dx 1 ^ ^ dxi 

where g = det(gjj) and (g tJ ) is the inverse matrix of (gij). We prove the following 

Theorem 4.1 Suppose the Ricci curvature Ric > —K for some K > 0, and suppose u is a positive 
solution of ( [^. j| ) with initial data uq > 0. Then 

2K 

\V\ogu\ 2 {t,x) < ^rlllognolU • (4.2) 

1 — e 2 
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The preceding theorem is proved by using similar computations as in the proof of Theorem 3.7 
but working on the orthonormal frame bundle O(M) over M. 

Recall that a point 7 = (x, e) G O(M), where (ei, • • • , e n ) is an orthonormal basis of the tangent 
space T X M at x £ M. Let 7r : 7 = (x, e) — > x be the natural projection from 0{M) to M. O(M) is a 
principal fibre bundle with its structure group 0(n). For the general facts on differential geometry, 
we refer to Kobayashi and Nomizu |15j . 

Suppose x = (x , • • • , x n ) is a local coordinate system on M, then it induces a local coordinate 
system 7 = (x fc , e*-) on 0(M) so that e 3 - = e*-^-. If L is a vector field, then L denotes the horizontal 
lifting of L to 0(M): 

L{x,e)=L\x)^-T{x%L\x)4-^ 

in a local coordinate system, where T^- are the Christoffel symbols associated with the Levi-Civita 
connection, and L = L l -J~. For a = 1, • • • , n and 7 = (x, e), then L a denotes the horizontal lifting 



dx 

of e a , that is 



The system {L\, ■ ■ ■ ,L n } is called the system of canonical horizontal vector fields. The mapping 
L : M n —> T(TO(M)) where = £, a L a , is defined globally, and is independent of the choice of a 
local coordinate system. Therefore 

n 

A 0(M) = y2 L a 
a=l 

is well defined sub-elliptic operator of second order on the frame bundle O(M), called the horizontal 
Laplacian. If / 6 C 2 then &o{M)f 71 = A/. For simplicity, any function / on M is lifted to a 
function / on O(M) defined by / = / o tt which is invariant under the group action by 0(d). 
The following relations will be used in what follows. 

L a /( 7 ) = Z a /o7r(7) = e ^, for 7 = (x^) GO(M). (4.3) 
We need the following geometric facts, whose proofs are elementary. 
Lemma 4.2 For a = 1, • • • , n we have 

-y n 1 " 

^0(M)L a - L a A ( M ) = - ^2[Lp, L a ]Lp + -'^2Lp[L/ 3 ,L a ]. (4.4) 

13=1 0=1 

Lemma 4.3 Suppose f is a smooth function on M and f = f o n is the horizontal lifting of f to 
O(M). 

1 ) For a, /3 = 1, • • • , n 

[L a ,L j8 ]/ = 0. (4.5) 

2) We have 

n 

{Lj){[L^L a ]Lpf) = Ric(Vf,Vf) (4.6) 

a,fi=l 
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where Ric is the Ricci curvature. 
3) We also have 

1 n 

&>o(M)L a f — L a Ao(M)f = 7;'^2[Lf3,L a ]Li3f. (4.7) 



2 

0=1 



Proof. The first identity (|4.5p follows from the torsion-free condition. To prove (j4.6j) we choose 

k 



a local coordinate which is orthonormal and dgij = (so that T^- = 0) at the point we evaluate 
tensors, so that 

I LLJ cVlc* 93 f 9T ^^ df 
t f t f - i * i d 3 / _ q i j'n% of 

P a Pl ~ 6a6 ^ dafldatdxi 6 P 6ae P dxi dx* ' 

and therefore 

which yields (|4.6p . (|4.7j) comes from (4.4) and (4.5). ■ 
Consider the following stochastic differential equation 



Mt)=Y t L a {i P (t))odw?, #)= 7 GO(M) (4.8) 



a=l 



on the classical Wiener space (fi, J 7 , P) of dimension n. The stochastic flow associated with (|4.8p is 
denoted by {y>(£, -,7) : i > 0}, which is a diffusion process in O(M) with the infinitesimal generator 

1 1 71 

0=1 

in the sense that 

ft 



M/ = f(t, ip(t)) - f(0, ip(0)) - J ^A 0(A/) /( S , <p( S ))ds 



2 

is a local martingale for any / 6 C 1,2 (IR+ x O(M)), and 

M / = E / (£«/)(*, • (4-9) 

We may express 

tp(t,w,j) = (X(t,w,-y),E(t,w,7)) 

where X(t,w,~f) = 7r((p(t, w, 7)), then {X(t, -,7) : t > 0} is a diffusion process on M starting from 
x = vr(7) with infinitesimal generator |A. {^(i, -,7) : i > 0} is a Brownian motion on M starting 
from x = 7r(7). 

In a local coordinate system (x fc ,e*), write 

¥»(t,-,7) = (^*(t,-,7);^(t,-,7)) 
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so that E a (t) = E l a (t)-^. Then the SDE (JOJ) may be written as 



dE)(t) = - Ea^8,fc=i rw^(t)£f(f) o d<. 



Let F(t) = (F(t)J) = E{t)- 1 . Then F(t)E(t) = I so that 



4.1 Gradient estimate 

Suppose now / satisfies the non-linear heat equation 

1 . <9\ . I 



(4.10) 



dFj{t)= Y, nXt) l 0J Fi(t)E^t)odwf. (4.11) 
a,P,l=l 

If / G C 1 - 2 (R+ x M), then 

/(T,X T ) = f(t,X t ) + j Yl E h a (s)^(s,X s )dw« 

+ [ T {lL + l A ) fis ' xs)ds ' (412) 

so according to (|4,3p 

f(T,X T ) = f(t,X t )+ [ T Y( L J)( s >( X s,E(s)))d< 

•** a=l 

+ r(£ + i A ) /(s ' x - )<is - (4 - i3> 



Let / be the horizontal lifting of / i.e. / = / ° vr, so that / satisfies the parabolic equation on 
0{M): 

(\^0 { M)-^f = -\Y\hf\ 2 . (4.15) 

v ' /3=1 

Let T > 0. Let = f(T - t, X t ) and 

Z? = YK(t)^L(t,X t ) = (Lj)(t,(X t ,E(t))) (4.16) 

for a = 1, • • • , n. Then 

n 

\Vf(T-t,-)\ 2 (X t ) = Y\Z?\ 2 (4-17) 



a=l 
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so that, according to (|4.12p . 

fT i rT n 



Y T = Y t + [ Z«dw a a -\ \ Y^\Z?\ 2 ds. (4.18) 



We will consider f|4. 18j) as a backward stochastic differential equation. 

Applying Ito's formula to L a f and the stochastic flow ip(t, -,7) to obtain 



dZ? = ^L f) (Lj)(t,<p(t,','y))dv>l 

(3=1 

d_ 



+ ( - A 0(M) + — ) LJ(t, <p(t, j))dt, (4.19) 



so that 



d\Z\ 2 = 2j2Ztd Z t+ E \LpLj\ 2 (t, V (t,;j))dt 

a=l a, /3=1 



n 



a=l 

However, by Lemma 4.3, 



2 J2 Z?{lplJ)dw P t + \LpLj\ 2 dt 

a,(3=l a, 13=1 

+2 Z? (~A (M) + §-) Ljdt . (4.20) 



1 d . , 

-A om + -\L 



(1 <9 \ 1 n 

(3=1 

n 1 n 



4 

,8=1 /3=1 



we therefore have 



n n 



d\Z\ 2 = 2^Z?^(LpLj)(dwe -Z?dt 

a=l (3=1 



2 

a, (3=1 a, (3=1 



+ £ \LpLj\ 2 dt + \ Y^{Lj){[L p ,L a )Lpf)dt. 

a,(3= 

By using Lemma 4.3 we obtain 

n 

d\Z\ 2 = 2J2Z?*£(L Lj)(dw?-zPdt') 

a=l (3=1 

n 1 
+ E \L p Lj\ 2 dt+ -mc(V f,Vf)(T-t,X t )dt. (4.21) 



n n 



a,(3=l 
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Define a probability Q by 



Rt where 



Rt = exp 



Then 



so that 



zW-i / \Z\ 2 ds 



Y T = Y t + f Z?dw« + \ Z t\ 2ds 
Jt *Jt a=l 



\ZJ ds 



Ft <4||/ | 



But on the other hand e 2 \Z t \ 2 is submartingale under Q, so that 



\ZJ 2 ds 



= J T ^{\Z s \ 2 \F t }ds 

Ks 

e 2 ds 

- x i z *i 

2 



which yields that 



and hence (14.211. 



1 — e 2 



K (T-t) 



00 



(4.22) 



5 Li-Yau's estimate 

In this section we prove Theorem 11.21 Thus, u is a positive solution to the heat equation 

1 d \ 

-A - — )u = on [0, 00) x M 

2 at/ 



(5.1) 



where M is a complete Riemannian manifold of dimension n, with non-negative Ricci curvature. 
Then / = log u is a solution to the semi-linear heat equation 



(5.2) 



Taking derivative in the parabolic equation (|5.2p with respect to the time parameter t one 
obtains 

& -Jt = Lf t (5.3) 



where 



df 



L = -A + V/.V . 
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By using the Bochner identity one can verify that 

d 



dt 



|V/| 2 = L|V/| 2 - iWjf - 2Ric(V/, V/) 



(5.4) 



Since A/ is the trace of the Hessian VV/ so that |VV/| 2 > ^ (A/) 2 , thus, since the Ricci 
curvature is non-negative, then 



Let 



||V/| 2 <L|V/| 2 -i(A/) 2 



G = -Af = \Vf\ 2 -2f t 



(5.5) 



By combining (|5.3p and (|5.5p together, we obtain 



0_ 

Of 



1 



G < LG - -G A . 



n 



(5.6) 



Next we apply the maximum principle to 

F -- 



*+i)(-A/) 



Since 



d_ 

at 



L F 



n C 



-G + 



n 

< i c 

n 



t 




1 




+ 




n 


C 


t 




1 




+ 




n 


c 







L G 



and therefore 



-G(l-F) , 

n 



By applying the maximum principle, we obtain the estimates in Theorem 11.21 

Corollary 5.1 Suppose u is a positive solution to 115.1]) , then the following Harnack inequality 
holds: 

u(t, x) 



u(t + s,y) 



< exp 



f t+s P ' 
/ K(a, ^)da 
Jt s 



where 



and 



Thus 



K(t,a)= sup \a^/4>(t,Y) -y\ 

Y:ip(t,Y)>0 J 



u(t, x) 
u(t + s,y) 



i;(t,Y)=2Y + j 
n 

/ 1 i 1 



c 



C + l 



< 



V 


~r(x,y) 2 ~ 


1 exp 


2s 





This follows by integrating the gradient estimates in Theorem 11.21 along geodesies, see [2] for 
details. 
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